Spurious Shell Closures in the Relativistic Mean Field Model 
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Following a systematic theoretical study of the ground-state properties of over 7000 nuclei from 
the proton drip line to the neutron drip line in the relativistic mean field model [Prog. Theor. 
Phys. 113 (2005) 785], which is in fair agreement with existing experimental data, we observe a few 
spurious shell closures, i.e. proton shell closures at Z = 58 and Z — 92. These spurious shell closures 
are found to persist in all the effective forces of the relativistic mean field model, e.g. TMA, NL3, 
PKDD and DD-ME2. 
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PACS numbers: 21.10.Dr, 21.10.Pc, 21.60.-n, 21.60.Jz 



The relativistic mean field (RMF) model has achieved 
great successes in studies of nuclear structure for many 
years. Recently, great efforts in the nuclear physics 

community have been taken to study exotic nuclei with 
extreme isospin ratios. In this respect, the RMF model 
is considered to be a very promising one due to its self- 
consistent description of the spin-orbit interaction. 0, 
A careful and systematic study of the ground-state prop- 
erties of over 7000 nuclei throughout the periodic table 
has been performed using the RMF+BCS model. Q 
It was shown that the RMF model can describe all the 
ground-state properties of those experimentally known 
nuclei very well. Nevertheless, a few deficiencies were 
observed: For example, the neutron skin thickness of 
208pj-, jg predicted to be 0.26 fm; while recent experimen- 
tal data seem to be favorable to a smaller value. This 
deficiency has been (partialljO removed by some lately 
developed effective forces. @, ^ Another deficiency, 
which has not received enough attention, is the appear- 
ance of spurious shell-closures at Z = 58 and Z ~ 92. 
In Ref. p|, nuclear binding energies around ^*°Ce and 
^^^U were found to be strongly overestimated. A closer 
examination of the corresponding single-particle spectra 
showed that they are caused by spurious shell closures at 
Z = 58 and Z = 92 (just beyond the conventional Z = 50 
and Z = 82 shell closures), which are further amplified 
by the co-occurrence of the = 82 and A^ = 126 neutron 
shell closures. 

In this letter, we investigate whether the same spuri- 
ous shell closures exist in other representative effective 
forces of the RMF model except TMA To save 

space, the formulation of the RMF model is not pre- 
sented here, which can be easilyfound in many review 
articles (see, for example, Ref. jj]). The effective forces 
we examine in this work are TMA, [HI which has been 
demonstrated to work very well throughout the peri- 
odic table, '6', "T] NL3,'T5| one of most-studied effective 
forces, PKDD, 9] and DD-ME2,_10J two lately developed 
density-dependent effective forces which have made some 
improvements compared to older ones. We believe that 
these four effective forces can represent most existing 
forces in the current RMF framework. Since an extensive 
study such as that performed in Ref. |^ is not possible for 



all the four forces, we would like to concentrate on -'^^"Ce 
and ^^^U, where the spurious shell closures at Z = 58 
and 92 were first clearly identified (see Fig. 2 of Ref. 0]). 
For comparison, we also study ^'^^Sn and ^"^Pb, which 
are used as reference nuclei. To avoid any possible am- 
biguity, all the results shown in this study are obtained 
with a spherical code and without the pairing correla- 
tion. This first approximation is due to the fact that 
these nuclei are found to be spherical,"^, and the sec- 
ond approximation (together with the first) implies that 
the so obtained binding energies for these nuclei are at 
the lower bound. 

In TableHl the binding energies of ^^^Sn, ^^^Ce, ^ospb 
and ^-'^^U obtained using the effective forces TMA, NL3, 
PKDD and DD-ME2 are tabulated. It is easily seen that 
all the four forces overestimate the binding energy 
by at least 8 MeV; while only TMA and PKDD overes- 
timate the binding energy of ^^°Ce by the same order. 
This is in agreement with the systematic study of Ref. 
6] . Among the four forces, the latest force DD-ME2 pro- 
vides the best description of the experimental data, but 
it still fails for ^^*U. This corresponds to the spurious 
shell closure at Z = 92, as we shall see in the follow- 
ing. It is interesting to note that PKDD does particu- 
larly well for ^'^^Sn and ^°®Pb while it is not so good for 
^^°Ce and ^^^U. Considering that these effective forces 
describe most nuclei throughout the periodic table with 



TABLE I: Binding energies of "^gn, ^osp^^ wo^g g^^^ 2is^ 
obtained with effective forces TMA, NL3, PKDD and DD- 
ME2 (in units of MeV). The experimental data are taken 
from Ref. [3]. 





«2Sn 


2°«Pb 


"OCe 


218u 


TMA 


1103.4 


1634.7 


1179.4 


1674.0 


NL3 


1103.2 


1637.5 


1176.3 


1674.7 


PK-DD 


1102.6 


1637.2 


1179.0 


1677.7 


DD-ME2 


1103.5 


1639.0 


1175.9 


1673.5 


Exp. 


1102.9 


1636.4 


1172.7 


1665.6 
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TABLE II: Proton shell closures &t Z = 50, 58 in "^Sn (""Ce) and those &t Z = 82, 92 in ^^^Pb (^^**U) obtained with effective 
forces TMA, NL3, PKDD and DD-ME2 (in units of MeV). The first column under ^^^Sn (""Ce) is the sheU closure at Z = 50 
and the second column is that at Z = 58. The first column under ^"^Pb (^^*U) is the shell closure at Z = 82 while the second 
column is that a.t Z = 92. 





^^^Sn 


^«'Ce 


20«Pb 


218u 


TMA 
NL3 
PKDD 
DD-ME2 


5.66 3.73 
6.15 2.81 
6.45 2.99 
6.45 2.39 


5.36 4.12 
5.78 3.50 
6.09 3.63 
6.01 3.08 


2.47 4.01 
3.50 3.20 
3.59 3.36 
4.06 2.74 


2.34 3.93 
3.14 3.54 
3.25 3.63 
3.64 3.28 


Exp. 


6.03 0.96 




4.20 0.90 





a discrepancy of 2-3 MeV. P. the 8 MeV discrepancy 
is rather large. This therefore indicates some physics 
that should not be ignored. Following Bertsch et al.flE^ 
these nuclei should be thought of as "critical nuclei" of 
the corresponding forces. 

In Ref. 6], the overbindings of ^^°Ce and ^^®U were 
attributed to the corresponding spurious shell closures 
at Z = 58 and Z = 92. In the following, we investigate 
whether it is still the case for NL3, PKDD and DD-ME2. 
The proton shell closures at Z = 50, 58 in "2Sn {^^°Ce) 
and those at Z = 82, 92 in ^ospb (2i8u) obtained with 
effective forces TMA, NL3, PKDD and DD-ME2 are tab- 
ulated in Table ^1 For comparison, the experimental 
data for ^^^Sn and ^°*Pb are also shown. For ^^^Sn, it 
is seen that NL3 agrees with experiment best; TMA un- 
derestimates the Z = 50 shell closure but overestimates 
the Z = 58 shell closure; while PKDD and DD-ME2 
are of similar quality. The large Z = 58 shell closure in 
TMA and PKDD persists even at ^'^'^Ce, which explains 
the overbinding of this nucleus as listed in Table For 
208pb, DD-ME2 agrees with experiment best; TMA un- 
derestimates the Z = 82 shell closure but overestimates 
the Z — 92 shell closure; while NL3 and PKDD are of 
similar quality. For ^^^U, all the effective forces predict a 
very large shell closure at Z = 92: it is even larger than 
that at 82 in TMA, NL3 and PKDD; while in DD-ME2 it 
is slightly smaller than that at 82. This extremely large 
Z — 92 shell closure is believed to be responsible for the 
overbinding of ^^^U as listed in Table It is interesting 
to note that the existence of these spurious shell closures 
seems to be a common feature of the current RMF forces 
(at least for those investigated here), while in nonrela- 
tivistic Hartree-Fock models, no such a common feature 
has been found. 

The nucleus is composed of both protons and neutrons. 
Therefore, one must study not only protons but also neu- 
trons to obtain a complete picture. In Table ITTll the neu- 
tron shell closures at TV = 82 in ^'^^Sn (^''"Ce) and those 
at = 126 in '^°^Ph (^^^U) obtained with effective forces 
TMA, NL3, PKDD and DD-ME2 are tabulated. It is 
seen that all the effective forces overestimate the A^ = 82 
shell closure by about 1.5 MeV. From ^^"^Sn to ^"^"Ce, 
this shell closure increases in NL3, PKDD and DD-ME2; 



TABLE III: Neutron sheU closures at iV = 82 in "^Sn ("°Ce) 
and those at = 126 in ^°*Pb (^^^'U) obtained with effective 
forces TMA, NL3, PKDD and DD-ME2 (in units of MeV). 





"^Sn 


""Ce 


208pb 


218U 


TMA 


6.78 


5.58 


3.58 


2.14 


NL3 


6.17 


6.38 


4.60 


3.07 


PKDD 


6.50 


6.65 


4.78 


3.16 


DD-ME2 


6.17 


7.07 


5.25 


3.72 


Exp. 


4.84 




3.43 





while it decreases in TMA. For ^ospb, TMA agrees with 
experiment very well; while all the others overestimate 
the experimental value with DD-ME2 having the largest 
discrepancy. From ^°^Pb to ^^^U, a decrease of about 
1.5 MeV is seen for all the effective forces . From Ta- 
bleland Table ITTll it can be concluded that the better 
description of the proton shell closure of •^"^Pb in NL3, 
PKDD and DD-ME2 compared to TMA is at the expense 
of the overestimation of the neutron shell closure, which 
is particularly true for DD-ME2. 

The importance of identifying these spurious shell clo- 
sures or "critical nuclei" can be summarized as follows: 
(i) It points out the regions in the periodic table that 
the present RMF model can not describe so well, (ii) It 
tells us to further study what causes these spurious shell 
closures and how to remove them. As a result, the cur- 
rent formulation of the RMF model can be improved. In 
the following, we discuss several approaches that might 
be helpful to remove the observed spurious shell closures. 
It should be noted that our discussions are restricted to 
the mean-field level and are not meant to be exhaustive. 
First, one may need a more balanced fitting strategy. 
For example, only the ground-state properties of twelve 
nuclei are fitted to obtain the latest effective force DD- 
ME2. Meanwhile these twelve nuclei are in some sense 
neutron rich. Since all the current effective forces in the 
RMF model adopt a similar procedure to fix their param- 
eter values, it is unsurprising that they all overestimate 
^^^U, which is proton rich. Although an extensive fitting 
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TABLE IV: Binding energies of ^^^^Sn, '^°^Ph, ^''°Ce and ^^^^U 
obtained with effective forces DD-ME2 and DD-ME2(M) (in 
units of MeV). The experimental data are taken from Ref. 





i^^Sn 


208pj^ 


"°Ce 


218u 


DD-ME2 
DD-ME2(M) 


1103.5 
1100.6 


1639.0 
1634.4 


1175.9 
1172.9 


1673.5 
1168.7 


Exp. 


1102.9 


1636.4 


1172.7 


1665.6 



strategy such as that used in obtaining the latest series 
of Hartree-Fock-Bogohubov mass tables 16] in the RMF 
model is not likely to be carried out in the near future, 
a more balanced fitting strategy can hopefully remove or 
at least reduce the overbinding of ^^^U and those around 
it. In this sense, taking into account the single-particle 
spectra, in particular those of ^°^Pb, during the fitting 
process might also improve the description of these crit- 
ical nuclei. 

Another alternative is to make magic nuclei slightly 
less bound. The reason is that particle-number conserved 
pairing methods usually give a few MeV of energy gain 
to even magic nuclei. |17l| Therefore, fitting the binding 
energies of magic nuclei exactly to the experimental data 
or even slightly larger than the experimental data (most 
current RMF forces behave in this way) leaves no room 
for further improvement and meanwhile overestimates 
magic number effects. Just to see how we can reduce 
the overbindings of ^""^Ce and ^^^U by making magic nu- 
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clei, particularly -'^^^Sn and ^''^Pb, less bound, we modify 
the effective force DD-ME2. It should be noted that the 
purpose is not to develop a new effective force, but to 
demonstrate our above statements. To keep things sim- 
ple, we only change the nucleon masses. The nucleon 
masses used in DD-ME2 are 939.0 MeV for both protons 
and neutrons. Til We reduce them to 938.5 MeV. The 
binding energies obtained with this modified DD-ME2 
(denoted as DD-ME2(M)) for ^^"^Sn, ^^^Ce, ^^^Pb and 
^^^U are tabulated in Table HVI One can easily see that 
now the discrepancies are more evenly distributed, i.e. 
serious discrepancies disappear. 

To make RMF calculations more reliable for astro- 
physical studies and/or studies of exotic nuclei, the pre- 
dictions of nuclear masses must be improved. It may 
or may not be possible within the current formulation 
of the RMF model. One may finally resort to a simi- 
lar procedure adopted in non-relativistic calculations.^^ 
Although recently constructed effective forces do show 
improvements compared to older ones, deficiencies re- 
main as demonstrated in this work. The recent system- 
atic study 6] showed that most discrepancies between 
the RMF results and existing experimental data actually 
originate from the regions surrounding the so-called crit- 
ical nuclei. Therefore, removing them should be viewed 
as an essential criteria for the next-generation effective 
forces. These works are underway. 
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